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Abstract 

The main purpose of the paper is an essentially probabilistic analysis of 
relativistic quantum mechanics. It is based on the assumption that whenever 
probability distributions arise, there exists a stochastic process that is either 
responsible for temporal evolution of a given measure or preserves the measure 
in the stationary case. Our departure point is the so-called Schrodinger prob- 
lem of probabilistic evolution, which provides for a unique Markov stochas- 
tic interpolation between any given pair of boundary probability densities 
for a process covering a fixed, finite duration of time, provided we have de- 
cided a priori what kind of primordial dynamical semigroup transition mecha- 
nism is involved. In the nonrelativistic theory, including quantum mechanics, 
Feyman-Kac-like kernels are the building blocks for suitable transition prob- 
ability densities of the process. In the standard "free" case (Feynman-Kac 
potential equal to zero) the familiar Wiener noise is recovered. 
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In the framework of the Schrodinger problem, the "free noise" can also be 
extended to any infinitely divisible probability law, as covered by the Levy- 
Khintchine formula. Since the relativistic Hamiltonians |V| and ^— A + m^ — 
m are known to generate such laws, we focus on them for the analysis of prob- 
abilistic phenomena, which are shown to be associated with the relativistic 
wave (D'Alembert) and matter-wave (Klein-Gordon) equations, respectively. 

We show that such stochastic processes exist and are spatial jump pro- 
cesses. In general, in the presence of external potentials, they do not share 
the Markov property, except for stationary situations. A concrete example of 
the pseudo differential Cauchy-Schrodinger evolution is analyzed in detail. The 
relativistic covariance of related wave equations is exploited to demonstrate 
how the associated stochastic jump processes comply with the principles of 
special relativity. 

1 The analytic continuation in time of holomor- 
phic semigroups as a mapping between two fam- 
ilies of stochastic processes 

1.1 Gaussian exercises 

The Schrodinger equation and the generahzed heat equation are connected by an- 
alytic continuation in time. The hnk (casually viewed as a kind of analogy or cor- 
respondence) can be implemented by a rotation in the complex time plane, taking 
the Feynman-Kac kernel into the Green function of the corresponding quantum me- 
chanical problem, which is an exploitation of properties of holomorphic semigroups 
generated by Laplacians and their sums with appropriate potentials. 

For V = V{x), X G -R, bounded from below, the generator H = —2mD'^A + is 
essentially selfadjoint on a natural dense subset of L^, and the kernel k{x,s,y,t) = 
[exp[—{t — s)H]]{x,y) of the related dynamical semigroup is strictly positive. The 
quantum unitary dynamics exp{—iHt) is a final result of the analytic continuation. 

As repeatedly emphasized @, 01) any temporal evolution that is analyzable in 
terms of a probability measure may be interpreted as a stochastic process. In view 
of the Born statistical interpretation postulate for quantum mechanics, the analytic 
continuation in time discussed above induces a class of probability measures, namely, 
consider p(x, t) = \ip{x, t)p as the density of a probability measure associated with a 
given solution ipi^x, t) of the Schrodinger equation. Then, it is possible to address the 
problem of that stochastic dynamics which would be either (i) measure preserving 
or (ii) induce the time evolution of the measure proper. Keep in mind that the 
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Schrodinger equation itself is not a genuine partial differential equation of probability 
theory; rather it is the Born postulate which embeds the unitary evolution problem 
into the probabilistic framework. 

A simple illustration of the analytic continuation in time is provided by consid- 
ering the force-free propagation, where the formal recipe gives rise to the equations 
of motion (one should be aware that to execute a mapping for concrete solutions, 
the proper adjustment of the time interval boundaries is indispensable): 



Then 



idtip = -DAtp — > dtd^ = DA9^ 

idt^ = DA^ — > dtO = -DA9, (1) 
it t. 

^jj{x,t) = [p^/^exp{iS)]{x,t) = I dx'G{x -x',t)-^{x',0), 



G{x - x\t) = {4mDt)-'/' exp[- ^^^-^], (2) 



9^{x,t) = / dx'k{x — x' ,t)6^:{x' ,0), 



k{x-x',t) = {ATrDty/'expl- ^\^], 
where the imaginary time substitution 

k{x — x' , it) = G{x — x' , t) 

k{x — x' , t) = G{x — x', —it) (3) 

seems to persuasively suggest the notion of "evolution in imaginary time", which in 
the usual interpretation relates quantum theory to an "imaginary time diffusion". 
Here we shall emphasize a different viewpoint in which the quantum dynamics and 
the so-called Euclidean dynamics (dealing with the Wiener process and condi- 
tional Brownian motions, for example) may both be seen as real time diffusions. 
At this point let us observe that given the initial data 

i;(x,0) = ina')-'/'expi-^) (4) 
the free Schrodinger evolution idtip = —DAip implies that 

^(x, t) = {-){a' + 2zDt)-'/' exp [ - (5) 
TT 2(a^ + ZtDt) 
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with 



p{y,0,x,t)p{y,0)dy (6) 



0, X, t) = (AnDtr'/^expl - ^^jjf^^^' ] 

where p{y, 0, x, t) is the (distorted Brownian) transition probabihty density for Nel- 
son' s diffusion ||. The transition density p{y,0,x,t) is not uniquely specified 
by (6), but the nonuniqueness problem for the diffusion process involved may be 
resolved when we consider transition densities for arbitrary intermediate times. 

For example 0, if for convenience we set = 2, Z) = 1, the transition density 
for two arbitrary times, t > s, reads 



p{y^ s, X, t) = [47r(t — s)] ' exp 



Ait-s) ^' 

c = c(.,t) = [ ^^-f t V , (7) 
1 + 

p{x,t) = / dyp{y,s,x,t)p{y,s), 



Notice that the s | limit of the transition density (7) does not coincide with the 
rescaled form of (6), |1 —t\ instead of (1 — t) appears in the exponent. It reflects the 
nonuniqueness of the definition of the transition probability density as long as we 
do not insist on having defined all intermediate densities as well 0. Anyway, while 
integrated with p{x, 0) they give the same output at time t > 0: 

2 2 

Po{x) = {27i)-'/^exp[~—] p{x, t) = [27r(l + t^)]-'/^exp[- ^ ^ '"^ 



Clearly p(x, t) admits a factorization p(a;, t) = \ilj{x,t)\'^ = 6(x, t)9*(x, t) . The stan- 
dard Madelung exponents R{x,t), S{x,t) such that ip{x,t) = exp[R{x,t) +iS{x,t)] 
are given by 

R{x, t) = -hn2nil + t') - j^^^, (9) 
Six,t) = ^ - ^arctan (t) , 
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and allow us to define the real functions 6(x, t) = exp[R{x,t) + S{x,t)], Q^:{x,t) = 
exp[R{x, t) — S{x, t)] which solve the pair of time adjoint generalized diffusion equa- 
tions 

dtO = -AG + ge, 

dtQ, = Ae, - QQ,, (10) 
ApV2 



Q(x,t) 



pi/2 



The diffusion governed by the pair of adjoint equations (10) belongs to the category 
of " Nelson's diffusions" |^ , but only in the framework presented here can it be singled 
out uniquely. It is exactly due to the "Schrodinger problem" uniqueness theorem 
0. It is really amazing that Schrodinger originated the problem of a stochastic 
interpolation between the prescribed input-output statistics data long before the 
modern probability theory was created. 

An interesting observation is that we can give the transition density (7) another 
form i, i 

p{y, s, X, t) = k{y, s, x, t) ^\ , (11) 

Q{y,s) 



lim 



AsiO 



^ [1 - y" k{y, s,x,s + /\s)dx] = Q{y, s) 



On the other hand, coming back to the previous notation and (5), we can straight- 
forwardly pass to 

e^ix,t) ■.= ^ix,-tt) = {-fl\a' + 2Dt)-"\xp[- ,f,„,J (12) 

TT Z\a + ZL>t) 

Let us confine t to the time interval [— T/2, T/2] with DT < a"^. Then we arrive at 

dt9, = DAe,, 

dtO = -DAO, 
T T 

where (we use an overbar to distinguish between the probability densities (14) and 
(6) or (7) respectively; notice also that 6 replaces 0) 
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with the interesting and certainly unexpected-if one follows traditional Brownian 
intuitions-outcome that: 

p{x,-t) = p{x,t) (15) 

for all \t\ < T/2. The density (14) refers to a conditional Brownian motion, and the 
interpolating probability density can be represented as the conditional probability 
density (identifiable as the Bernstein density P]) 

-( ,N _ -"o,a;,t)A;(x,t,0,ao) 



/c(0, -ao,0,ao) 

'4:D{t-s) 



k{y,s,x,t) := [4nD{t-s)]-'/'exp[-^^—^], (16) 



«2 

" 2D 

Since p(x, t) trivially factorizes into a product of solutions of time adjoint heat 
equations (set 9{x,t) = [k{0, —ao,0, ao)]~^^'^k{x,t,0, ao)), we have in hand the mi- 
croscopic transport recipe 

- 6(x t) 

p{y, s, X, t) = k{y, s, x, t) ' (17) 

with the heat kernel kiy, s, x, t), (16), and the property (to be compared with (11)): 
HmAsj,o-^-[l — J k{y, s,x,s + As)dx] = (18) 

The resemblance of formulas (11) and (17), (18) is not accidental, and suggests 
that any given Feynman-Kac kernel can be used to generate a probability measure. 
In addition, we should keep in mind that the two levels of probabilistic description 
i.e. (11) and (17), are indirectly linked by the analytic continuation in time of a 
holomorphic semigroup with the Laplacian as its generator. 

Remark 1: The time developement of a density p{x,t) (p respectively) is dic- 
tated by the Fokker-Planck (second Kolmogorov in the probabilistic lore) equation. 
It is instructive to notice that idfi/j = —Aijj upon setting v = 2Re^ and u = 2Im^ 
gives rise to dtp = — V(fp), which may be rewritten as dtp = Ap — V(6p) with 
b = u + v. Proceeding analogously with p = k{xi,ti, x, t)k{x, t, X2, t2)/k{xi,ti,X2, ^2) 
where k{y,s,x,t) is the heat kernel, and s < t, ti < t < t2, while ti,Xi,t2,X2 are 
fixed, we immediately recover dtp = Ap — V(6p) with b = b{x, t) = 2Vk{x, t, X2, ^2)/ 

k{x,t,X2,t2). 
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Furthermore to this remark, let us emphasize that the emergence of the nonvan- 
ishing drift field b{x, t) is not connected with any external force . It is a traditional 
assumption when studying the Brownian motion in a conservative force field to 
define the drift as being proportional to the force itself (Stokes law); these Smolu- 
chowski diffusions form a subclass of problems we are considering P|, however the 
previous drift is an exclusive effect of the conditioning. 

1.2 The Schrodinger problem: from Feynman-Kac kernels 
to probability measures 

The previous examples are very particular solutions of what we call Q the Schro- 
dinger problem of deducing the probabilistic interpolation (stochastic process) con- 
sistent with a given pair of boundary measure data at fixed initial and terminal 
time instants ti <t2. Originated by Schrodinger himself 0], the problem was solved 
much later ^ by invoking the machinery of Bernstein stochastic processes; see 
also Ref. 4. For our purposes the relevant information is that [0], if the interpolating 
process is to display the Markov property, then it has to be specified by the joint 
probability measure [A and B are Borel sets in R): 



JR 

where we assign densities to all measures to be dealt with, and the density m{x, y) 
is given in the functional form 



involving two unknown functions f{x) and g{y) which are of the same sign and 
nonzero, while k{x,s,y,t) is any bounded strictly positive (dynamical semigroup) 
kernel defined for all times ti < s < t < t2. The integral equations (19) determine 
functions f{x),g{y) uniquely (up to constant factors) in this case 0. 



By denoting 6^{x^t) = J f{z)k{ti,z,x,t)dz and 9{x,t) = J k{x,t, z,t2)g{z)dz it 
follows i, S I, i that 






(19) 




m{x,y) = f{x)k{x,ti,y,t2)g{y) 



(20) 



p(x, t) = 9{x, t)9^{x, t)= p{y, s, X, t)p{y, s)dy, 



(21) 
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p{y, s, X, t) = k{y, s, x, t)^f^, 

h < s <t <t2 

Hence the transition densities are intrinsically entangled with the dynamical semi- 
group kernels in the solution of the Schrodinger stochastic interpolation problem. 
The crucial step in the construction of any explicit propagation consistent with the 
boundary measure data is to decide what is the appropriate dynamical semigroup. 

We shall address the issue in its full generality. Strictly positive semigroup 
kernels generated by Laplacians plus suitable potentials are very special examples 
in a surprisingly rich encompassing family. First of all, the concept of the "free 
noise" , normally characterized by a Gaussian probability distribution appropriate to 
a Wiener process, can be extended to all infinitely divisible probability distributions 
via the Levy-Khintchine formula ^. It expands our framework from continuous 
diffusion processes to jump or combined diffusion-jump propagation scenarios. All 
such (Levy) processes are associated with strictly positive dynamical semigroup 
kernels, and all of them give rise to Markov solutions of the Schrodinger stochastic 
interpolation problem (19)-(21). 

Remark 2: Apart from the wealth of physical phenomena described in terms 
of Gaussian stochastic processes, there is a number of physical problems where the 
Gaussian tool-box proves to be insufficient to provide satisfactory probabilistic ex- 
planations. Non-Gaussian Levy processes naturally appear in the study of transient 
random walks when long-tailed distributions arise [1T0| , ITTj |T2l . They are also found 



necessary to analyze fractal random walks [|T3], intermittency phenomena, anoma- 
lous diffusions, and turbulence at high Reynolds numbers |jlO|, |15|, 

Let us consider Hamiltonians of the form H = F{p), where p = —iV stands for 
the momentum operator and for — oo < k < +00, F = F{k) is a real valued, bounded 
from below, locally integrable function. Then, exp{~tH) = exp[—tF(k)]dE(k), 
t > 0, where dE{k) is the spectral measure of p. 

Most of our discussion will pertain to processes in one spatial dimension, and let 
us specialize the issue accordingly. Because {E{k)f){x) = J^oo^^Pi''-P^)fip)^Pj 
where / is the Fourier transform of /, we learn that 

/+00 
exp{-tF{k))dE{k)f]{x) = 
-00 

I (^xp[~'tF{k)] — [ exp{ipx) f {p)dp\dk = (22) 

J — 00 CirC J — oo 



1 r+oo 



exp{—tF{k))exp{ikx)f{k)dk = [exp{—tF{p))f{p)Y{x) 

00 
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where the superscript V denotes the inverse Fourier transform. 

Let us set kt = -^[exp{—tF{p)]'^, then the action of exp{—tH) can be given in 
terms of a convolution: exp{—tH)f = f * kt, where (/ * g){x) := J^g{x — z)f{z)dz. 

We shall restrict consideration only to those F{p) which give rise to positivity 
preserving semigroups: if F{p) satisfies the celebrated Levy-Khintchine formula, 
then kt is a positive measure for alH > 0. The most general case refers to a con- 
tribution from three types of processes: deterministic, Gaussian, and an exclusively 
jump process. We shall concentrate on the integral part of the Levy-Khintchine 
formula, which is responsible for arbitrary stochastic jump features: 

/+00 j^py 
[exp{ipy) - 1 - — — ^]iy{dy) (23) 
-oo L -\- y 

where z/((iy) stands for the so-called Levy measure p, [T7[] . 

The disregarded Gaussian contribution would read F{p) = cf. Refs. 3-7 

for an exhaustive discussion of related topics. In this case we know in detail how 
the analytic continuation in time of the Laplacian generated holomorpic semigroup 
induces a mapping to a quantum mechanical (since the Schrodinger equation is 
involved) diffusion processes. 

Our further attention will focus on two selected choices for the characteristic 
exponent F{p), namely: Fq{p) = \p\ and Fm{p) = y/p'^ + m'^ — m,m > 0, where we 
have chosen suitable units so as to eliminate inessential parameters. (The relativistic 
Hamiltonian is better known in the form -y/m^c^ -|- c^p^ _ -^j^gj^g q jg the velocity 
of light.) 

The respective Hamiltonians (semigroup generators) Hq, are pseudodifferen- 
tial operators. The semigroup kernels kt, k"^ in view of the "free noise" restriction 
(no potentials, will be defined in below) are transition densities of the jump (Levy) 
processes regulated by the corresponding Levy measures u^^dy), Umidy). The af- 
filiated Markov processes solving the Schrodinger problem (19)-(21) immediately 
follow. It is instructive to notice that as in the case of Gaussian derivations (1), it 
is the case 6{x,t) = 1, 9^{x,t) := 'p{x,t) for which the pseudodifferential analog 
of the Fokker-Planck equation, as a consequence of [exp{—tH)'p\{x) = p{x,t) and in 
view of the identification F{p — iV) := H takes the fundamental form 

Fo{p) ^ dtp{x,t) = -\V\p{x,t) (24) 

or 

FM =^ dtpix, t) = -l^-A + m^ - m]p{x, t) (25) 

respectively. Let us emphasize that the existence and uniqueness of solutions proof 
for the Schrodinger problem extends to all cases governed by the infinitely divisible 
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probability laws, and has nothing to do with the "nonrelativistic" or "relativistic" 
options. The particular choice of semigroup generators, which are called "relativistic 
Hamiltonians" links the standard Schrodinger problem discussion with relativistic 
dynamics. But only after an analytic continuation in time, unless only stationary 
problems are studied (see the forthcoming discussion). 

Although the pseudodifferential generator of the semigroup implies that the 
Fokker-Planck equation is no longer exclusively differential but an integro-differential 
equation, each solution p{x, t) in the present case is nevertheless a solution of a 
partial differential equation of higher order. Specifically, the respective partial dif- 
ferential equations are of the second order: 

Fo{p) =^ Oep{x, t) = (At + A)p(x, t) = 0. (26) 

Alternatively, if we set p{x,t) := p{x,t) exp{mt) in (25) then: 

Fmip) =^ (At + A)pix,t) = m'pix^t) 

^ (27) 
(-□^ + m2)p(x,t) =0 
where dtp = —\/—A + m'^p holds true instead of (25). 



Our two semigroups are holomorphic [|TB| , hence we can replace the time parame- 



ter t by a complex one a = t + is, t > so that exp{—(jH) = exp{—(jF{k)) dE{k). 
Its action is defined by 

[exp{-aH)]f = iifexpi-aF)]'' = f*K (28) 

to be compared with (22). Here, the kernel reads k„ = --^[exp{—aF)Y. Since H is 
selfadjoint, the limit t I leaves us with the unitary group exp{—isH), acting in the 
same way: [exp{—isH)]f = [fexp{—isF)Y, except that now his := -^[exp{—isF)Y 
in general is not a measure. In view of unitarity, the unit ball in is an invariant 
of the dynamics. Hence density measures can be associated with solutions of the 
Schrodinger pseudodiferential equations: 

Fo{p) =^ iM{x, t) = I V|V^(x, t) (29) 

or 



F^{p) =^ idtipix, t) = [^J-A + m2 - m]ilj{x, t) (30) 
provided with the appropriate initial data functions ^/'(x,0). 
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An obvious consequence of (29), (30) is that the partial differential equation of 
the second order (26) takes on a familiar relativistic form 

Fo{p) =^ a^ix, t) := (-A + At)^{x, t)=0 (31) 

while after setting ip{x,t) = ipixjt) exp{imt), we arrive at the Klein-Gordon equa- 
tion: 

Fmip) ^ (□ + m2)^(x, t) = (32) 

where the D'Alembert operator □ = — A + replaces its Euclidean counterpart 
-De in (27). 

We have thus reached a point, at which the main questions addressed in the 
present paper can be precisely stated: 

(i) What are the stochastic processes consistent with the probability measure 
dynamics p{x,t) = |?/'(x,t)p determined by pseudodifferential equations (29) and 
(30)? 

(ii) Can we extend the Schrodinger problem idea to the special relativistic domain 
and be able to reproduce the interpolating stochastic process from the given input 
p{x,ti) and output p{x,t2), ti < t2 statistics data, just as in the nonrelativistic 
(Laplacian generated motion) case ? 

(iii) To what extent can we attribute a definite probabilistic meaning to solutions 
of the relativistic wave equations (31), (32) ? 



2 Can we associate Feynman-Kac kernels with 
the pseudodifferential Schrodinger dynamics ? 

Given the Schrodinger equations (29), (30). To set them in the Schrodinger problem 
framework of Section 1.2 we need to choose any normalized solution and then take 
the associated probability density p(x, t) := \ipix, t)p as the boundary data at times 
ti < ^2- However, as stated before some additional requirements must be met, 
specifically the Markov property is necessary |^. One should keep in mind that if 
we do not insist on the Markov property for the interpolating process, then a solution 
of the problem involves the general Bernstein processes [0 for which a reformulation 
in terms of a pair of time adjoint generalized diffusion equations no longer exists. 

We have chosen two rather special pseudodifferential counterparts of the Lapla- 
cian guided by two reasons: (i) their similarity on analytic grounds (the same criteria 
p!8| for the existence of the bound state spectrum if summed with suitable poten- 
tials, which we shall need in the sequel), (ii) the claim of Ref. |]19[ that the pertinent 



stochastic process in the mass m > case actually displays the Markov property. 
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If the Markov property would hold true for the relativistic Hamiltonian generated 
dynamics, we would be able to repeat almost all steps of the previous Schrodinger 
problem analysis |, However, the situation is not that simple, and the 

subsequent argument excludes the Markov property, in all nonstationary situations, 
in a flat contradiction with general statements by De Angelis . 



Before embarking on this issue, let us introduce some probabilistic notions, which 
will tell us how to work with pseudodifferential operators. We shall notice that for 
explicit computational purposes, the Cauchy generator | V| is much more suited than 
the m > relativistic Hamiltonian. It is a real disadvantage when dealing with Levy 
processes that rather limited number of concrete examples is available, in contrast 
to the wealth of the general theory. 

The Levy-Khintchine formula (23) tells us that the action of the Hamiltonian 



{Hf)ix) = - lj,f{x + y)-f{x)-'-^^]v{dy) (33) 



H = F(— iV) on a function in its domain can be represented as follows p, [17 

yVfix) 

It is important to observe that for the "free noise" processes whose semigroup 
generators are |V| and v^— A + — m we do know explicitly their kernels (tran- 
sition probability densities) and the involved Levy measures, as well as about the 
extension of the Feyman-Kac path integral construction of the semigroup kernels 
to these particular Levy processes |jl8|, |21|, in case of arbitrary space dimen- 
sions. Therefore we feel free to use the Feynman-Kac kernel notion instead of the 
semigroup kernel. 

For the Cauchy process, whose generator is |V|, we deal with a probabilistic 
classics [g 0: 

'^ = h ^ ^ ^' = l it-syl'lx-yy ^ ^'^^ 

<s <t 

{exp[ipX (t)]) := / exp{ipx)p{x,t) dx = exp[—tFo{p)] = exp{—\p\t) 
The characteristic function of k^{y, s, x, t) for y, s fixed, reads exp[ipy— \p\(t—s)], and 



the Levy measure needed to evaluate the Levy-Khintchine integral reads p3|, ^ ^ : 

■h%0,0,y,t)]dy=^ 
t ny 



uo{dy) := Umtio[-k° {0,0, y,t)]dy = — (35) 



In the case of the relativistic generator -^Z— A + — m, formulas determining 
the stochastic jump process are much less appealing 



{exp[ipX (t)]) := exp[—tFm{p)] = exp[—t{J p"^ + w? — m)] 
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m texp{mt) 



p{x,t) = --^^^^ K,{mV^^) (36) 

[exp{-{t - s)Fm{-iV))]{x -y) = k"'{y,s,x,t) := p{x -y,t- s) 

in 

Um{dy) = ——Ki{m\y\)dy 

where Ki{z) is the modified Bessel function of the third kind of order 1. 

We are interested in acting with the pseudodifferential generators H = F(— zV) 
on functions in the exponential form (recall the familiar Madelung procedure in the 
Gaussian case) f{x,t) = exp^{x,t): 

^\/ \ /" r \ \ y^' {x)exp^{x) 
{Hexp(pj{xj = — J [exp'P{x + yj — exp9[x) — ^— j- — {I'ldy) = 

r y(i>'(x) 
= exp^{x) / \exp{^{x + y) - ^{x)) - 1 - - — ^-^]iy{dy) (37) 

JR 1 + y 

where <l>'(x) = V$(x). Since (if$)(x) = - Jj^[^{x+y)-^{x)-y^' (x) / {l+y'^)]iy{dy), 
we can make a safe rearrangement of (37): 

(i7exp$)(x) = exp<!>{x)[{H^){x) - f {exp^^y - 1 -^.^y)v{dy)] (38) 

:= ^{x + y) -^{x) 

In application to the pseudodifferential dynamics idtip^Xjt) = {Hip){x,t) with 
ip = exp{R + iS), we shall investigate its implications for the real functions G = 
exp{R + S) and G* = exp{R — S); our argument will admit a trivial extension from 
H to H + V situations. 

Remark 3: Experience p^ , ^ with the Gaussian (standard Laplacian gener- 



ated) noise proves that the Madelung substitution ip{x^t) = exp[R{x,t) +iS{x,t)] 
would associate with the Schrodinger equation a pair of time adjoint generalised 
diffusion equations where the Feynman-Kac potential (time dependent in the gen- 
eral case) equals ^;^[2Q{x,t) — V{x)]. Here Q{x,t) = 2mD'^^j^(x,t) and V{x) 
is taken as an external conservative force potential. Let us emphasize that V{x) 
actually was the Feynman-Kac potential of the dynamical semigroup prior to the 
analytic continuation in time procedure. The mapping V{x) — >• 2Q{x,t) — V{x) is 
an effect of the analytic continuation in time, as manifested on the level of the as- 
sociated Feynman-Kac kernels. Previous comments suggest that the nonrelativistic 
formalism can be viewed as a kind of probabilistic reinterpretation of Bohm's point 



of view Specifically, the function —Q(x,t) is the familiar de Broglie-Bohm 
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'quantum potential". The analogous connection is generally invalid in the context 



of the Klein-Gordon equation, as explained in Ref . [|3^ . 

To this remark, let us add that the analytic continuation in time as outlined 
in Section 1 is not a well known procedure although repeatedly mentioned in the 
previous publications It is new as a regular method. Also, it is not a part of 

Nelson's "quantum fluctuations" point of view, just as the Feynman-Kac kernels 
were not ingredients of Nelson's theory. The idea is new, though developed on 
the basis of Zambrini's [Q] and Carmona's ^ research. Refs.2-6 give a complete 
description of the state of art in this respect. Therefore we do not propose to indicate 
any further connections with Nelson's stochastic mechanics than by referring to 
Nelson's monograph . The paper itself proposes a new, self-contained probabilistic 
analysis which, in addition to Zambrini's work is motivated by the paper due to De 
Angelis . 



In view of (38), the pseudodifferential Schrodinger equation idtipix, t) = Hip{x^ t) 
implies the following time evolution of the Madelung exponents: 

dtR = HS - / [exp{Ra:y) sinS^y - S^y]diy{y) 

J R 

dtS = -HR + / [exp{Rxy) cos Sxy - 1 - Rxy]diy{y) (39) 

J R 

where H = F[—i'V). 

By employing (38) with respect to p^^^ = exp{R), we arrive at: 



Q '■= = HR- I [exp{R,y) - 1 - R,y]du{y) (40) 

^ R 



Hp'/' 
and hence: 

dtS = -Q+ f exp{Rxy) [cos{S,y) - l]diy{y) (41) 

R 

The same procedure can be repeated for 6 = exp{R + S) and 9* = exp{R — S), 
where equations (39), (41) imply: 

dt& = HQ + 6 [-2(5 + / exp{Rxy)[-sin S^y + cos S^y + expiS^y) - 2]du{y)] 

<J R 

(42) 

dtO* = -HQ^ + Q^[2Q - / exp{Rxy)[sinSxy + cos S^y + exp^-S^y) - 2]du{y)] 

J R 

In contrast to the Gaussian case [^, |^, equations (42) do not take the form of a 
time adjoint pair, unless some additional restrictions are imposed on the Madelung 
exponent S{x^ t) (notice that we have restored time dependence, skipped before 
for convenience). An obvious demand is S{x + y,t) = S{x,t) for all y,t, and any 
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fixed X. But then, equations (42) would manifestly refer to the stationary (measure 
preserving) random dynamics, governed by the pair of equations: 



dtO = He- 2QQ 

dtQ, = -HQ, + 2QQ, (43) 

which are mutually time adjoint. Hence they would fall into the Schrodinger prob- 
lem framework, with a trivial implication that the measure preserving process is 
Markovian. This however cannot be a property of the "free" dynamics since we 
need external potentials to secure stationarity. Let us therefore make an essential 
amelioration by performing the previous analysis for the case idtip = [H + V)il! with 
V = V{x). Then, the stationary system of equations (43) would take the form: 

dt& = HQ- {2Q + V)Q (44) 

dS, = -HQ, + {2Q + V)e, 

which upon substituting S{x,t) = Et, where E is a. constant, yields a pseudodiffer- 
ential version of the Sturm-Liouville problem: 

Hp^/\x) - [2^^ + V{x) - E]p'/^{x) = (45) 



pl/2 



V{x) - E = - 



p'/^{x) 



to be solved (for a chosen value of E) with respect to the square root of the proba- 
bility density p(x), once the external force potential V{x) is selected. 

This problem has its Gaussian counterpart in the study of the measure preserving 
dynamics 2^ , and in the present context it can be solved by invoking those 



potentials for the original pseudodifferential Schrodinger equation, for which the 
bound states (i.e., stationary solutions) have granted the existence status. The 
relevant analysis has been carried out in the studies of the relativistic stability of 
matter ||26|, |18], |2^, In addition we know [^, |18|] that in the stationary case, the 
Feynman-Kac path integral generalization to Levy semigroup kernels is available. 

However, the Markov property cannot automatically be attributed to the non- 
stationary dynamics, as described by (42). Below, we shall make a careful analysis 
of the Cauchy - Schrodinger {H = |V|) dynamics, to produce a definite counterex- 
ample, for which the unrestricted equations (42) would hold true, but the associated 
random dynamics would be non-Markovian. 
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3 The Cauchy-Schrodinger dynamics 



3.1 Construction of an explicit nonstationary solution 

While it is clear that exp{—t\V\) and exp{—it\V\) have a common, identity operator 
limit as i J, 0, an analytic continuation of the Cauchy kernel by means of (28) gives 
rise to: 

1 1 is 

kis = ::[Six -s) + S{x + s)] + -V- 



2^ ^ ' ' ^ ' ' 'k' — 
Here, we use the usual notation for the Dirac delta functionals, and the new time 
label s is a remnant of the limiting procedure t | in o" = t+is. The function denoted 
by is j 'nix?' — s^) comes from the inverse Fourier transform of —-^^sgn{p)sin[sp). 
Because of 

[sgn{p)Y = i\l^V{-) (47) 

V TT X 

where V{-) stands for the functional defined in terms of a principal value of the 
integral. Using the notation 5±s for the Dirac delta functional 5{x^ s): 

[sin{sp)Y = i^iS, - S_,) (48) 

we realize that 

-V^ = f (49) 

TT — ZTT X 

is given in terms of the implicit convolution of two generalized functions. 

Let iIj{x,0) = f{x) := be a normed function, which we take as the 

initial data for the Cauchy-Schrodinger evolution. 

With the unitary kernel kis in hand, we can define the pertinent evolution in 
terms of a convolution ijj{x, s) := f * kis = * (^s — ^-s) * ^(^)- Let us consider 



* /.i,,P(i) = ^pr/(i^i^rf,= (50) 



/ 



27r a; 27r 

i , . , (-^ f(x — s — v) , f°° fix — s — y) , , 
—lim,^o[ -dy+ ^ ^dy 

ZTT J-oo y Je y 

Because of: 

^1 ^ — j^iij [ , =1 — A(x — s) arctan (x — s — y) 

(51) 
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where 

^ = 

l + (x-s 



^ = TTVZ—X^^ (52) 



the definite integrals in (50) read: 



J-oo V V TT 



y V TT ^1 + - s + ey 



71" 

{x — s) arctan {x — s + e) + —{x — s)] (53) 
fix — s — y) , /2 , . , e 



Je V V TT 



and therefore: 



y V TT ^1 + (x - s - e)2 

TT 

— {x — s) + {x — s) arctan {x — s — e)] 
p /•+~ /(a; - s - y) _ V^{x - s) 



^-oo y 1 + [x - sy 

By proceeding analogously (with replacing 5s in (50)) we find: 

W-oo 1/ " 1 + (X + 5)2 ^^^^ 

All that finally implies (remembering that f{x) — /(I + x^) — ip{x, 0)): 
^{x, s) = [expi-isH) f]{x) = ^ [/(x + s) + /(x - s)] + 

[{x - s)f{x -s) - {x + s)f{x + s)] (56) 
with an interesting formula for the time developement of the probability density: 

p{x,s):= |V^(x,s)|^ = {1 + s"^) Jpo{x + s)po{x - s) (57) 



poix) = mx,0)\' = [f{x)]' = ^ 



TT (1+^2)2 

Notice that by a direct evaluation, we can check the normalization identity p{x, s) 
dx = 1. 

Now, we can address the problem of whether the stochastic process, implying 
the propagation (57) of the probability density, is a Markov process. 
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3.2 The nonstationary Cauchy-Schrodinger stochastic pro- 
cess is not Markov 

3.2.1 Candidates for the transition probability density 

We are interested in representing the time evolution of p{x,t), (57), in the integral 
(probabilistic transport rule-looking) form: 

,(.,t)^ f p(y,0,.,t)Mv)dy^ f P(y,s,.,t),(y.s)dy (58) 



without bothering at the moment whether we can assign to p{y, 0, x, t) or p{y, s, x, t), 
s < t, any true meaning of the transition probability density of a stochastic process. 

Since p{x, t) is a probability density, we can evaluate its characteristic func- 
tion (Fourier transform |^): 4>{p,t) := ^/27ip{p,t). Of course, (f>{p,0) = f{p) = 
/^[l/7r(l + x'^)]exp{ipx) dx = exp{—\p\) is a characteristic function as well. 

By observing that {■ip{x,0) := ipoix)) 

\ + i) + - *)i = /f (,'11)^/4 (59) 

while 

we arrive at (cf. (56) for the definition of tlj{x,t)) 



In view of (61), we have 



p{p, t) = — / exp{ipx) p{x, t) dx 



27T Jr 

^ exp (ipx) [Reip{x,t) — - Imip{x,t)]dx (62) 



27r JR t 



where 



So that 



ip{x,t) = — / exp(±iqx — \q\ — i\q\t) dq (63) 
\/27r JR 



^(j){p,t) = p{p,t) = exp {-\p\) [cos {t\p\) + ^sin(t\p\)] (64) 



271 V27r t 
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Finally, in view of po(p) 



(l/\/2^)(l + \p\)exp{-\p\), we find {t > 0): 



Po{p) 



cos{t\p\) + ^sin(t\p\) 

1 + IpI 



(65) 



and {0 < s <t): 



pip, t) 



p{p, s) 



cos{f\p\) + j-'^in{f\p\) 
cos{s\p\) + ^sin{s\p\) 



(66) 



Ignoring the issues of existence and positive-definiteness of Fourier transformed 
integrands, we can proceed in the standard way: 



where the formal, homogeneous in space "transition probability density" p{y, s, x, t), 
s < t, trivially satisfies the formal Chapman- Kolmogorov identity: p{y, s, x, t) 
— /i?P(?/) u)p{u, z, X, t)dz , s < u < t. 
To go beyond formal arguments: 

(i) We need to prove that the function p{y,0,x,t), t > is a well defined tran- 
sition probability density of the stochastic process transporting the initial (time 0) 
density into the terminal (time t >0) one. 

(ii) We need to demonstrate that the would-be transition density p{y, s,x,t), s < 
t, is a well defined probability measure, and actually we shall prove that it is not , 
which excludes the Markov property for the stochastic process under consideration 
in agreement with our previous conclusions of Section 2. 

3.2.2 Existence of the probabilistic transport from time to time t >0 

In agreement with (65)-(67), we can introduce an integral kernel p(|/, 0, x, t) effecting 
the transport of p{y, 0) into p{x, t) , t >0: 



At the moment, its status as the transition probability density of a stochastic process 
is not settled. We must know whether its Fourier transformed integrand is positive- 
definite function and satisfies a normalization identity JjiP{y, 0,x,t) dx = 1 for all 
times t > 0. 





(67) 




(68) 
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2 

.sin{tp) y 1 fW 



In view of the homogeneity in space we observe that p{y, 0, x, t) = p{0, 0,x — y, t) 
and so we pass to the notation 

V2^pix, t) = + ^£!!^]v (69) 

The function 

9ix) ■= -4=( i ^ I y ix) = TT / exp(-ipx) (70) 

V27r 1 + IpI 277 Ji?^ 1 + IpI 

will play a distinguished role in what follows. Indeed, because 

[cositp)r = ^l^{6t + 5-,) (71) 

tp ' tV2^[-*'*l 
where = 1 if < t, and = if |x| > t, we can rewrite p{x,t) as follows: 

p{x,t) = ^[g{x + t) + g{x-t)]+ (72) 

In view of (72), properties of p{x,t) are completely determined by those of g{x). 

Eq.(70) suggests that g{x) itself might be a probability density. For this to hold 
true, 1/(1 + \p\) must be a characteristic function, and we are in the framework 
covered by the classic Bochner Theorem P, ^ Our integrand 1/(1 + \p\) is 
continuous on R and equals 1 when \p\ = 0. It is well known ^ that such a function 
is a characteristic function if and only if it is positive definite. To be a positive 
definite function h{p) must satisfy the inequality 

^lJ,=lh{p^ -Pj)XAj > (73) 

for any finite sequence of complex numbers Ai, A2, A„, any sequence of points 
Pi,P2, ■■■,Pn in R , and any n = 1,2, .... 

The identity (73) is trivially satisfied by h{p) = 1/(1+ |p|), and as a consequence it 
is a characteristic function of the probability density g{x), (70). Then Jj^g{x)dx = 1 
follows. 

The function p = p{x, t) is real and even in x, which allows us to consider a 
nonnegative semiaxis x > for the moment. Let < x <t. Because Jj^g{x)dx = 1, 
we find that ^ ^ 

g{y)dy < 1 ^ - - -(g * X[-t,t]){x) > (74) 
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which imphes that 



p{x,t) = ^[g{x + t)+gix-t)] + ^ - l-{g*x[~t,t]){x) > (75) 



If X > t it suffices to consider 



1 g{x + t) + gix - t) 

g{y)dy < 



2t Jx-t 2 

^[g{x + t) +gix-t)] - ^^{g * Xht,t]i^) > ^ pix,t) >0 (76) 

Hence, p{x,t), (72) is strictly positive on R, and moreover it has a unit normal- 
ization. This implies that p{0, 0,x — y, t) is a well defined probability density. 

Let us add that lim\x\^ooP{x,t) = and li'mtioP{x,t) = 6{x). At the points 
X = ±t, p{x,t) develops singularities in view of lim^f^o g{x) = oo. The location of 
the singularities is quite significant, since they make a clear distiction between the 
random propagation regimes with jumps of size less than t, and those of size greater 
than t, for each terminal time instant of the evolution of p{x,t). 

Notice also that as a suitable probability density p{x, t) leads to several finite 
moments: /fj xp(x, t)dx = , x'^p{x, t)dx = t^. On the other hand, the long-tailed 
Cauchy transition density (34) has no finite moments at all. 

3.2.3 Violation of Markov property 

In view of our previous discussion we have specified consistent distribution functions 
for the process, which in principle permits its construction ( chap. 15.). For the 



Markov process what needs to be specified [|, ^ |32| are the transition probabihties 
for all intermediate times of the considered evolution and the initial distribution 
(which is given in our case). We shall investigate the existence of the consistent 
transition probability densities for the process in question. 

If the Cauchy-Schrodinger process is to be Markov, the integral kernel: 

( ^\ ^ [ \- ( M cos(t|p|) + \sm{t\p\) 

p{y,s,x,t) = — / exp[tp{y-x}\ . , I . . , dp (77) 

zn JR cos{s\p\) + ^sm[s\p\) 

must be the Fourier inversion of a characteristic function. A necessary condition for 

cos{t\p\) + \sin{t\p\ 

Kp, s, t) = . I |. 1 . . I |. (78) 

cos[s\p\) + -sm[s\p\) 

to be a characteristic function is the positive-definiteness of h{p,s,t) as a function 
of p E R for all intermediate time instants < s < t. 
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We shall prove that for each terminal time instant t we can single out earlier time 
instants s, such that h{p, s, t) is not a positive-definite function. As a consequence, 
the intermediate propagation cannot be Markov. 

To this end, let us notice that for a fixed t, and s > we can rearrange the 
denominator in (78): 

cos(s|p|) + -sin(s\p\) — cos(s|p| — as) (79) 

s cos as 



tana 



1 



s 



Notice that 



, , 1, 2N + 1 ^ , , , 

IpI = -{as H tt) =^ cos{s\p\ - q;^) = (80) 

S Zi 

for all integer iV. 

Our h{p,s,t) should satisfy the condition (73). Let us choose the simplest case 
of n = 2 in this formula, and consider llfj^'Ji[pi — pj)XiXj > 0. Because of (78), the 
two by two matrix hij :— h{pi — pj) , = 1, 2 has matrix elements hn = 1 = /i22 
and hi2 — /i2i = M with: 

I .cos{t\pi-p2\) + \sin{t\pi-p2\) 

M = cos{as) ^ 7 (81) 

cos[s\pi -P2\ - Ois) 

To have det[hij\ > we need \M\ < 1. This condition can always be violated 
by choosing any pair pi,p2 for which, at a given time s, the numerical value of 
\pi — P2\ — \p\ is close to any of those introduced in the formula (80). 

Therefore, the considered stochastic process (with the transition mechanism (67)) 
is not Markov, as anticipated on the basis of arguments of Section 2. 



4 Meaning of the pseudodifferential stochastic prop- 
agation: an insight into jump features of the 
process 

4.1 Fokker-Planck equations 

The probability density p{x,t) (respectively p{x,t)) was a fundamental entity in all 
our previous considerations: either (i) providing the input-output statistics data 
for the Schrodinger-random dynamics reconstruction problem, or (ii) providing the 
time evolution of the probability measure for the whole time interval of interest, so 
that the transition probability densities could be sought for. 
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In the Gaussian case of Remark 1 we dealt with the temporal evolution of 
the probability density given in its traditional Fokker-Planck form appropriate for 
Markov diffusion processes. In connection with the pseudodifferential ("free noise") 
dynamics, equations (24), (25) are an obvious extension of the previous notion to 
a class of jump processes. We shall extend the usage of the name Fokker-Planck 
equation to any first order in time differential equation determining the space-time 
properties of p{x,t) or p{x,t). 

Let us investigate the time developement of p(x, t) = 9{x, t)9^{x, t), where 9{x, t), 
6^:{x, t) come out as solutions of the temporally adjoint pair of equations of the form 

dtO = He-Ve (82) 

with the initial/terminal data f{x), g{y) of the Schrodinger problem (19), (20) and 
a Feynman-Kac potential V. Then, in view of (37) and 6 = exp{R + 5), 6^ = 
exp{R — S), we get an evolution equation for the probability density: 

dtp{x, t) = 9, {x, t) {HO) (x, t) - 9{x, t) {Hd,) {x,t) = (83) 

/ [-9,{x, t)9{x + y,t) + 9{x, t)9,{x + y,t) + 2p{x, t)VS{x, t)—^]du{y) 
jR i- + y 

Following the traditional recipes when dealing with Levy measures [|], let us consider 
an open neighborhood of the origin |e| ^ 1. Instead of integrating over all possible 
jump sizes, let us integrate over jumps of size \y\ > e > 0). The removal of this 
lower bound as e will eventually amount to evaluating the principal value of 
the integral. In case e > 0, we can safely remove the compensating term including 
y/ {l+y"^) from the integral, and restrict considerations to the contribution from the 
first two terms only. 

Our purpose is to establish a connection between (83) and the conventional 
theory of jump stochastic processes, as developed in Integrating over a Borel 
set A C R, X E A we get: 

[ dx [ [-9^{x,t)9{x + y,t) +9{x,t)9^{x + y,t)]di^{y) = 

JA J\y\>e 

f , [ / X r _/ ^(^ix + y^t) _ 9{x) 

dx XA{x)[-p{x,t) + p{x + y,t)-——-\du{y) = 84 

JR "'|?/|>e 9[x,t) 9{x + y) 

f dxp{x,t) f ^(^ + t) ^^^^-^ + y) - Xa{x)] du{y) 
JR J\y\><i o[x,t) 
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where we interchanged the order of integrations, and made appropriate adjustments 
of integration variables {x ^ x — y and y —y), while exploiting the property 
du^—y) = —du{y) of measures (35), (36); Xa{.x) is an indicator function of the Borel 
set A (Z R, equal to 1 when x & A and otherwise. 

In the present case, (82), we deal with a Markov process with transition probabil- 
ity densities given for arbitrary time instants: p{x, t) = J^piy, s, x, t)p{y, s)dy, s < t. 
By invoking the standard wisdom about jump Markov processes and exploiting 
limtisP{y, s, A,t) = XA^y), for any Borel set A C R away from (— e, +e), we can 
define the jump process running with jumps of size \y\ > e > 0. It should be viewed 
as an approximation of the original stochastic process governed by (83), with the 
initial data p(x, 0) common for both: 



dtMAt) = / qix,t,A)p^{x,t)dx + {v)Ait) / -J^duiy) (85) 

Jr J\y\>el + y^ 

where ^ 

q{x,t,A) := limuit:^^-—^[p{x,t, A, u) - Xa{x)] = 

^(^ + ^'^) [^^[x + y)-XA{x)]du(^y), (86) 

\y\>e 0{X,l) 

{v)A{t) := / p{x,t)[2VS{x,t)]dx 



Here q{x,t,A) > for all x which are not in A, in agreement with |32|. We have 
also introduced a pseudodifferential counterpart of the current velocity field v{x, t) = 
2VS{x,t), previously attributed to diffusion processes (cf. Remark 1), where the 
probability conservation law (a continuity equation in another lore) dtp = — V(fp) 
plays the role of the Fokker-Planck equation. 

Notice, that in the particular case of 6{x, t) = 1 for all x, t, and = 0, Eq.(82) re- 
duces to the "free noise" situation covered by the Fokker-Planck equations (24), (25). 
Then, q(t,x,A) = J\y\>e[XA{x + y) - XA{x)]diy(y) , while R = -S ,p = exp{2R) = 9^ 
implies {v)A(t) = — p(x, t)|^ where [a,b] := A C R. 

Now, let us address the Fokker-Planck equation for the pseudodifferential- Schro- 
dinger dynamics case, which we consider in the form analogous to (82); see also (1) 
for comparison: 

idt^ = Htlj + V'ip (87) 
idtip = —Hip — Vip 

We re-emphasize that to define the probability density p{x,t) = [^/^(x, t)p one actu- 
ally employs solutions of the time adjoint pair of Schrodinger equations. 
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In view of (87), the pseudodifferential continuity equation follows: 

dtp{x,t) = -i[^p{x,t){Hip){x,t) - ^{x,t){H^){x,t)] = (88) 

-i I [-i^{x,t)ij{x + y,t) + ij{x,t)i^{x + y,t) + 2ip{x,t)VS{x,t)-^]du{y) 
Our next step is a repetition of the procedures behind (84), which implies: 
dtp{x,t)= f [2I[,lj{x,t)^{x + y,t)] + 2p{x,t)VS{x,t)-^]du{y) 

JR i- + y 

dx / 2I['ip{x,t)'ijj{x + y,t)] = 
dx [ XAix)2p^/^{x,t)p^^'^{x + y,t)sin[S{x,t) - S{x + y,t)]diy{y) = (89) 

R J\y\>e 

f p^^'^ix + ?/) 

p{x,t)dx / sin[S{x + y,t) - S{x,t)][xA{x + y) -XA{x)]diy{y) 

R J\y\>e pV^(x) 

where 2[f{x,t)] stands for an imaginary part of a complex function f{x,t). So, a 
counterpart of (85) reads: 

dtp,{A,t) = [ q{x,t,A)p,{x,t)dx+{v)A{t) f -^du{y) (90) 

Jr J\y\>e 1 + r 

where, however 

q{x, t, A) := f J[^1^±M] [xj,{x + y)- XA{x)]dv{y) (91) 

J\y\>e t{x,t) 

no longer can be derived from transition probability densities of the process, as in 
the previous discussion (86), because in general our process is not Markovian. At 
least in the case of nonstationary dynamics, the only transition probability density 
which is at our disposal connects an initial instant of the evolution with any later 
one. In fact, we might even not be sure that q{x, t, A) is a well defined probabilistic 
object, because of the presence of sin[S{x + y,t) — S{x, t)] in the integrand. At this 
point an observation of [|1^ helps. Namely, in view of the identity: 



dx \ip{x + y,t)^lj{x,t)\[xAix + y)-XAix)]du{y) = (92) 

R J\y\>e 

valid for Borel sets A (Z R, which are away from (— e, +e), we can always pass from 
(89) to the rearranged form of (91): 

implying that q{x, t, A) is positive for all x which are not in A, as should be the case 
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4.2 The jump processes toolbox 

To have a better insight into stochastic jump processes associated with the Fokker- 
Planck evolutions (83), (85) and (88), (90) respectively, some further knowledge of 
the general theory is necessary. We shall try to minimize the level of sophistication 
by invoking arguments based on the exploitation of the standard Poisson process. 

A random variable X taking discrete values 0, y, 2y, 3y, with y > is said to 
have Poisson distribution V{X), A > with jump size y, if the probability oi X = ky 
is given by P{X = ky) = -^exp{—X). The characteristic function of V{X) reads: 

\k 

E[exp{ipX)] = exp[X{e'Py - 1)] = S^e"^ e'^^^^^ = j:^e"'^Py^ P{X = ky) (94) 

and its first moment equals E[X] = X. Notice that P{X = 0) = exp{—X), hence 
the numerical value of A > tells us what is the probability of a jump not to occur 
at all for a given Poisson process. For the Poisson random variable with values 
b + ky, k = 0,1, ...we would get 

E[exp{ipX)] = exp[ibp + A(e*P^ - 1)] . (95) 

If we consider n independent random variables Xj,l < j < n such that Xj has Pois- 
son distribution V{Xj) with jump size yj, then a new process X can be introduced 
with the distribution of Xi + ... + X„ so that its characteristic function reads 

E[exp{ipX)] = exp[S^=iAj(e^P» - 1)] (96) 

The exponent in (96) would include an additional term ipJ^^bj if nonrandom shifts 
of each jump kyj by bj would be allowed. 

We can admit not only jumps of fixed magnitudes yi, ...,yn but also jumps cov- 
ering an arbitrary range in R^. Let the distribution function of the magnitude of 
the jump be P{x < y) = fi{y). A possible generalization of (96) to this case is 

E[exp{ipX)] = exp[ f (e^^^ - l)dii{y)\ (97) 

assuming that the integral in the exponent exists. Notice that (96) is recovered, if 
we set 

d^^{y) = Tr^^^X,5{y-y,)dy (98) 

The convergence of the exponent in (97) may be jeopardized in cases when jumps 
of very small amplitude are allowed to occur very often, while we take for granted 
that jumps of very large size seldom happen. On the other hand for any Borel 
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set A C R bounded away from the origin, the process Xa of jumps bounded by A, 
has the characteristic exponent /^(e*^^ — l)dii{y), and the expected number 
of jumps of size A is equal to Ijl{A). We can say that the processes of jumps of 
different sizes proceed independently of one another, and the jump process of jumps 
of size [y,y + Ay), Ay ^ 1 contributes a Poisson component with exponent function 
approximately equal to (e*^^ — l)fj.{[y,y + Ay)). At the moment the processes have 
only upward jumps, hence their sample paths are nondecreasing. 

For a process with the characteristic exponent — (23), we can consider its 
restriction to upward jumps of size y > e > 

<Pt{p) = / [e'^' - 1 - T^]diy{y) = / [e'^' - l]du{y) - tpbt (99) 
•'y>e i-+y •'y>e 

Clearly, we deal here with a process of the type considered before, and might try to 
isolate contributions from jumps of size [y, y + Ay) by considering a coarse-graining 
of a Borel set A of interest. A formal substitution of (98) in (99), with d/j, replacing 
di/, gives rise to 

E[exp{ipX)] = expm^^[Xj{e'Py^ - 1) - ip -^^] (100) 

to be compared with (95). 

Further specializing the problem to relativistic Hamiltonians, we notice that 
the corresponding Levy measures uo{y) and u^niy), (35), (36) are even under space 
reflections, hence dv{—y) = —dv{y) in these cases. Consequently, we can easily 
extend our discussion to jumps of all sorts in R, i.e. y can take values in both i?+ 
and with the only restriction to be observed that |y| > e > 0. Notice that we 
shall deal with two processes, which run separately with either positive or negative 
jumps, and there is no common jump point for them. This fact means that they are 
independent components of the more general process defined by: 

Mp)=f [e'"" -l]dv{y)-iph (101) 

where in the case of the Levy measures given in (35), (36) the deterministic term 
identically vanishes in view of 

= ht + K = [ j^Mv) + / Tx;;^^^(^) ^ (102) 

Jy>e l + y^ Jy<-e I + y^ 
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All our steps (94)-(101) involved the fact that we deal with infinitely divisi- 
ble probability laws. One additional important property about them is that if 
ea;p0(p) is a characteristic function of a given probability distribution, then [ea;p0(p)]* 
= exp[t(j){p)], t > is likewise a characteristic function of an infinitely divisible 
probability law again. This feature readily extends our discussion to time-dependent 
stochastic processes (time homogeneous with independent increments, associated by 
us with the "free noise"). Obviously, for such processes E[exp{ipX (t))] = exp[t(f){p)] 
while EA[X{t)] = tui^A), and our previous arguments retain their validity with re- 
spect to 

E[exp{ipX{t))], = exp[t(p,{p)] = exp[t [ (e^^^ - l)du{y)] (103) 

J\y\>^ 

dtp{x,t) = -{Hp){x,t) =^ dtp^{A,t) = dx[ [xA{x+y)-XA{x)]du{y)]p^{x,t) 

JR J\y\>e 

Remark 4: In fact, (102) means that the Fokker-Planck equations (85), (90), 
if specialized to Levy measures (35), (36), involve exclusively the integral term on 
their right-hand-side: 

dtMAt)= I q{x,t,A)p^{x,t)dx (104) 

JR 

dtp^{A,t) = / q{x,t,A)p^{x,t)dx 

JR 

where an overbar distinguishes between probabilistic quantities characterising differ- 
ent families of stochastic jump processes (86) and (91) respectively. Let us emphasize 
that the simplification (104) occurs only in the |?/| > e > jumping size regime. The 
real role of two terms in (102) is to compensate the divergent contributions from 
the Levy measure when the principal value integral e — limit is considered; then 
the standard jump process theory (104) does not apply. Anyway, those two terms 
are irrelevant for any e > 0, irrespectively of how small e is. 

One might expect that an infinitesimal (jump size) surrounding of the origin 
gives a dominant jump contribution to the process. However, generally it is not the 
case: explicit solutions (34), (36) and (57), (72) indicate that for times t > the 
leading contribution does not necessarily come from jumps of infinitesimal sizes. 

5 Relativistic wave equations and associated sto- 
chastic processes 

We mentioned before that solutions of our pseudodifferential-Schrodinger equations 
solve the relativistic wave (or matter-wave in the Klein-Gordon case) equations as 
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well; see (29)- (32). Since each particular solution has an undoubted probabilistic 
significance, we can re-analyze the old-fashioned problem ^3, 34 1 of a" single-particle 



interpretation" for free Klein-Gordon solutions and analyze the same problem for 
the D'Alembert equation solutions, from a novel perspective; see also [|^. As 
well, we can benefit from relativistic covariance properties of wave equations to 
understand how the pseudodifferential-Schrodinger stochastic processes comply with 
the principles of special relativity. 

To begin with, let us consider the Klein-Gordon equation for a particle of mass 
m > 0: 

(□ + m2)0(f,t) = (105) 

The spacetime metric signature is diag{g^^) = (1, — 1, — 1, — 1, ), and the system 
of units is ^ = c = 1. In view of the polar (Madelung) decomposition of the 
complex wavefunction, (f){x,t) = exp[R{x,t) +iS{x,t)], we can split (105) into two 
real equations: 

{d,S)id>'S)=m^ + ^^ (106) 
d,f = 

r := ^^im - <pd^^] = -Pid'-S) 

where p{x,t) = |0(x, t)p = exp[2R{x,t)]. 

We can handle the m = case corresponding to the D'Alembert equation in the 
same way, and the only change in formulas (106) would be the absence of the 
contribution. 

We have noticed before, (32), that if ipi^x, t) is a solution of the pseudodifferential- 
Schrodinger equation idtip = [^/—A + m'^ — m]tp, then ilj{x,t) = il){x,t)exp{—imt) 
is a positive energy solution of the free Klein-Gordon equation (□ + m?)ip{x,t) = 
0, since we surely have idtip = a/— A -|- m'^tp. It is clear that the time adjoint 
Schrodinger equation refers to negative energy solutions of the Klein-Gordon equa- 
tion. Notice that we need both positive and negative energy solutions to create 
(upon normalization) a probability density p(x, t) = ip{x, t, )ip{x, t). 

Remark 5: At this point it is useful to emphasize that the time-like component 
j'^{x,t) of the current j'^^x^t) is not a probability density itself; by wrongly |3^, |3^ 



and per force assuming that it generally would be the case, all known paradoxes and 
difficulties underlying the refutation of the Klein-Gordon equation as the proper rela- 
tivistic generalization of its nonrelativistic Schrodinger partner are revealed [^, ^ . 
The positive energy spectrum is not just correlated with positive (negative) val- 
ues of jo{x,t), although one can establish such a correlation for the total "charge" 
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e / jo{x,t)(fx |3^, ^ by assuming that ejo{x,t) is interpreted as the "charge" den- 
sity. Even then, a clean partition of the positive and negative energy spectra into 
sets associated respectively with particles and antiparticles distinguished by the sign 
of the charge density is impossible. 

Remark 6: The subject of our considerations is essentially a probabilistic anal- 
ysis of relativistic quantum mechanics In view of our previous discussion it is 



clear that p{x, t) = \tp{x, t)p is a probability density of a well defined stochastic pro- 
cess, which is non-Markovian in nonstationary situations. Then, for a general Borel 
set A (Z we have defined a measure p{A, t) telling us what is the probability for a 
jump to have its size matching a point y G A, at time t. We are not inclined to think 
that a concrete jump refers to an actual "physical particle" that jumps in space. In 
nonrelativistic quantum mechanics, a standard interpretation of p{x,t){Ax)^ as a 
probability to locate a physical particle in a cube of volume (Ax)^ seems to be con- 
sistent. On the contrary, in relativistic quantum mechanics, the notions of position 
and localizability and their relation to any experimental determination of the phys- 
ical particle position have been a subject of vigorous disagreements and no general 
consensus has been reached. An acceptance of the Newton- Wigner localization in 
the configuration-space approach to relativistic quantum theory implies the general 
breakdown of causality 0, ^ ^ and inevitably implies superluminal 



effects (instantaneous spreading of a localized wave packet). It is by no means a sur- 
prise if one carefully looks into the jump process features as revealed in the present 
paper. Let us also mention that the Newton- Wigner localization of mass m = par- 
ticles is in general impossible (well known exceptions are massless particles of spin 
and massless spin 1/2 particles possessing two helicity states) , while we know how 
to assign the probability density notion, and hence a probability measure p{A, t) to 
a class of solutions of the D'Alembert equation, see e.g. the arguments of Section 3, 
albeit possibly with no connection to any ^ "position operator" notion. 

Each scalar positive energy solution 0(x, t) of the free Klein-Gordon equation 
(105) can be represented in the manifestly Lorentz covariant form: 



4>{S, t) = J d^k e^-*^"^") 5(A;^F - m')e{ko) $(A;o, k) (107) 

where k := {ko,k), k^k^ := fcg — fc^, $(/c) is a scalar and Q{ko) is the Heaviside 
function equal to IJf ko > and to O otherwise. The representation (107) extends to 
all solutions of idtij = y/—A + m'^ip, and upon changing ko — > —ko in 6(/co) followed 
by a complex conjugation of (108), to solutions of the time adjoint equation as well. 
It implies that general solutions of those pseudodifferential-Schrodinger equations 
form Lorentz invariant subspaces in the linear space of all solutions to the free 
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Klein-Gordon equation. 

However, we can not directly infer from the above facts any information about 
how a given pseudodifferential-Schrodinger stochastic jump process is perceived by 
different relativistic observers. For example the normalization of Schrodinger wave 
functions is not a relativistically covariant notion. 

At this point we adopt the standard definition of the Klein-Gordon scalar product 



d^x [0ia/-A + m202 + (^-A + m20i)02] (108) 

which is independent of the specific space-like surface of integration. Both positive 
and negative energy solutions are covered in this definition, albeit separately, with 
no superposition. The integrand in (108) should be compared with the time-like 
component j^{x) of the conserved four-current j^{x), Eq.(106). 
Since the Newton- Wigner position operator 

— * 

[fw0](x) = z[Vg - -J^—mk,, k) (109) 



ko = y k'^ + rn^i j = 1, 2, 3, see (107), is Hermitean with respect to the scalar prod- 
uct (108), one can introduce a covariant localization notion for all positive en- 
ergy solutions of the free Klein-Gordon (and hence pseudodifferential-Schrodinger) 
equation ^3|, |19|. Indeed, given a positive energy solution of the free Klein-Gordon 



equation 0(a;), which we know to solve the pseudodifferential equation as well, then, 
we can introduce a new solution for both of those equations as follows 



0(x) -> (t)Nw{x) ■■= [(-A + m2)V40](a;) (no) 

If we take for granted the Klein-Gordon scalar product normalization, Eq.(108),of 
0(x), we realize that the common solution of the Klein- Gordon and pseudodifferential- 
Schrodinger equations: ip{x, t) := 4>nw{x, t), may be consistently normalized accord- 
ing to the Schrodinger equation rule: p(x) := |'?/'(a;)p =^ Jj^^s d^xp{x,t) = 1. As a 
consequence, and in part because of this normalization, we have succeeded to asso- 
ciate the previously investigated stochastic jump process with the Newton-Wigner 
localization. Clearly, we deal with the probability measure identifying the probabil- 
ity that the Newton-Wigner "particle" can be found at time t in the spatial volume 
A, with the probability of spatial jumps bounded by this volume, at time t: 

Prob[X{t) e A] = f p{x,t)d^x= [ \[{-A + m^Y/^(l)]{x,t)\^d^x (111) 



A 



The inhomogeneous orthochronous Lorentz mapping x' = Ax + a {x'^ = K'^x^ + 
a'^, Aq > 1) can be associated with the scalar transformation rule for Klein-Gordon 
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wave functions 0'(xg,x') = 0(xo,x) =^ 4''{x) = (j){A~^{x — a)). The transformation 
acts invariantly in positive and negative energy subspaces of solutions, respectively. 
This fact implies an extension to pseudodifferential-Schrodinger equations of motion. 

The Klein-Gordon equation is form invariant: (□' + m^)0'(x') = (□ + m^)0(x) = 
0, which allows us to associate with (j)'{x'), while normalized according to (108), a 
pseudodifferential-Schrodinger stochastic process according to (110): 

i^'ix') := [(-A' + m2)0f /4(x') 



idt'^\x,t') = + ni^i>'{x\t') 

^ (112) 

idt'^'{x',t') = [^-A' + m2 -m]^'(f ,t') 

ip'{x,t') := exp{imt')ip'{x ,t') 

The transfer of data about the pseudodifferential-Schrodinger process from one 
inertial observer to another is completely determined by Eq.(112). 

Now, let us consider the first of equations (87) for the choice H = y/—A + 
of the Hamiltonian: [idt — (V — m)]ilj = \/—A + m'^ip, where we can regard a 
general potential V — m as a time-like component of a four-vector. In case of 
electromagnetic interactions, the presence of the (— m) term can even be attributed 
to a gauge transformation — > + d^x with x = —mt and A = 0. In particular, 
the relativistic stability of matter studies of the existence of bound states for the 



pseudodifferential-Schrodinger Hamiltonians, involve the Coulomb static potential 
eAo{x) = — (Ze^/r), r := (We recall that pseudodifferential Hamiltonian 

spectral problems are not limited to the electrostatic potential only, but their range 
of applicability is much wider |[T8[] .) 

Complex conjugation converts the forward equation into its time adjoint and it 
is clear that [~idt ~ {V — m)] ip = y/—A + ip holds true. However, no immediate 
connection with the general form of the Klein-Gordon equation in the presence 
of electromagnetic interactions (charge e particles), {id^ — eAn){id'^ — eA'^)(f){x) = 
m?(j){x)^ can be established, in general. 

On the other hand, a pedestrian intuition behind the associated notion of a 
relativistic atom is quite helpful for a deeper understanding of the particular role of 
the Lorentz covariance of the Klein-Gordon and D'Alembert equations in the context 
of free pseudodifferential Schrodinger equations. Namely, the atom itself is always 
considered to be at rest, with the frame of reference attached to the nucleus, which 
in turn is a source of an electrostatic field. In this particular frame of reference the 
pseudodifferential-Schrodinger equation [idt — {V — m)]ilj = a/— A + m'^ip is defined. 
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The same pertains to the general pseudodifferential-Schrodinger problem with a 
minimal electromagnetic coupling due to external fields 



[idt - (ev4° - m)]tp{x, t) = ^S3^J(zVj - eAj){^^ - eA^)] + m? ^(x, t) (113) 

which is a frame-of-reference-dependent notion, see e.g. ||48|| . 

Remark 7: A Euclidean version of (113) was investigated in |^ and an ex- 
plicit construction was given of the Feynman-Kac path integral formula for the 
corresponding semigroup kernel. It involves paths, and conditional measures over 
paths, of a time homogeneous Levy process. 

Remark 8: The problem of how a stochastic process can be perceived by differ- 
ent relativistic observers has been considered before ||4^ in connection with certain 
(Markov) rotational diffusions on an 5*3 manifold {SU{2) x SU{2) case specialized 
to spin 1/2), with the Euler angles parametrization established relative to a fixed 
three-dimensional orthonormal basis. 

Our discussion was confined to the mass m > case, but in view of the general 
non-existence of the Newton- Wigner localization in the mass m = case, an ex- 
tension of previous relativistic covariance arguments needs some care. The present 
localization is known to be admissible for massless spin particles. As well, we 
do not literally need the Newton- Wigner-like position operator and the associated 
notion of localization at a spatial point to invoke a substantial part of the previous 
arguments. (In fact, a maximal localizability of photons on a circle, hence in a 
subset of R^, was established in [0.) As long as (108) is extended to mass m = 
particles as the normalization condition for wave functions, and the Borel set A is 
never point-wise, we can safely go through (110)-(112). 

Remark 9: There have been numerous attempts to associate the Klein-Gordon 
equation with stochastic processes. In addition to |1^ let us mention a number of 



50| , p^ , p3| , |5J, |5^, ^ ^ |36|. None of them can be viewed as 



other attempts |£l], 

a "derivation" of the Klein-Gordon field from certain "first" stochastic principles. 



Except for [|T^, all these attempts exploited a formal similarity of Eqs.(106) to local 
conservation laws shared by nonrelativistic Markov diffusions and to analogous 
laws in relativistic kinetic theory The status of the Markov property 



has been found disputable, since its violation is implicit if the relativistic invariance 
of diffusion (Kolmogorov) equations in Minkowski space is required 
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In the present paper, we have found the Markov property admissible only in case of 
the measure preserving (stationary case) stochastic jump dynamics; see e.g. Section 
2. 

Acknowledgement: We express our gratitude to Professor G. G. Emch for a 
helpful discussion. 
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